Apéndice A 


Nociones básicas sobre Cálculo Tensorial 


A.1. Introducción 


El concepto de tensor tiene su origen en la evolución de la geometría diferencial de Gauss, Rie- 
mann y Christoffel. La necesidad del cálculo tensorial, conocido también como cálculo diferencial 
absoluto, como rama sistemática de la matemática, se debe a Ricci y a su discípulo Levi-Civita, que 
publicaron en colaboración el primer trabajo sobre esta materia: Métodos del cálculo diferencial ab- 
soluto y sus aplicaciones, en "Mathematische Annalen", vol.54 (1901). 

El objeto principal del cálculo tensorial es la investigación de las relaciones que permanecen 
invariantes cuando se cambia de un sistema de coordenadas a otro. Las leyes de la física no pueden 
depender del sistema de referencia que elija el físico con fines descriptivos. Por eso, es estéticamente 
deseable y muchas veces conveniente, utilizar el cálculo tensorial como fundamento matemático en 
que se puedan formular tales leyes. Einstein, en particular, lo consideró un excelente instrumento para 
la presentación de su teoría general de la relatividad. 

El cálculo tensorial alcanzó gran importancia y es hoy día indispensable en sus aplicaciones en la 
mayoría de las ramas de la física teórica. 


A.2. El convenio de sumación de Einstein 


Una suma cuyos sumandos se obtengan dando los valores 1,2,...,n a ciertos subíndices de su 
término general, se indica comúnmente con el símbolo Y junto con la indicación del intervalo de 
variación y la forma del término general. Por ejemplo: 


n 
aibí +42b3+-->+4pb, = Y aibi 
i=l 
n 
a1bjci +azbjc2 ... Anb ¡Cn = Y aibjci 


i=l 
En el cálculo tensorial aparecen con mucha frecuencia sumas del tipo de los dos ejemplos anteriores 
en que la suma se verifica respecto de dos subíndices repetidos de su término general. En este caso es 
cómodo y abrevia mucho la escritura convenir que, para tales subíndices, se suprimirá el símbolo de 
la suma, estableciendo así el siguiente convenio: 
Cuando en una expresión monomia figuren dos subíndices repetidos, se entenderá que se trata de 
una suma en la que los subíndices repetidos van sumados de 1 a n. 
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La letra n hace referencia a la dimensión del espacio y su valor resulta siempre claro del tema gue 
se está tratando. Es conveniente ser cuidadoso porgue cuando no se guiera expresar la suma habrá gue 
indicarlo explícitamente. 

Usando el convenio de Einstein los ejemplos anteriores se escribirían como: 


n 
Y aibi = aib; 
i—I 


n 
y aib jci = aib jc; 
i=1 


A.3. Representación de un vector en un sistema cartesiano 


En primer lugar veamos cuál es el sistema de referencia que vamos a emplear para representar a 
los vectores. Vamos a distinguir entre dos tipos de sistemas de referencia atendiendo a la orientación 
relativa de sus vectores básicos: 


= Directo: se llama así al S.R. que tiene sus vectores básicos orientados de tal forma que el pro- 
ducto mixto! de los tres resulta ser 1: 


ölni 


Dicho de otra manera, cuando un sacacorchos girando del primer vector al segundo por el 
menor ángulo avanza en la dirección y sentido del tercer vector. Este tipo de sistema también es 
Ilamado en ocasiones diestro o dextrógiro. 


= Inverso: Se trata del caso opuesto al anterior. Ahora tenemos: 


e (ene) =-1 


o si un sacacorchos girando del primer vector al segundo por el menor ángulo avanza en la 
misma dirección pero sentido opuesto al del tercer vector. Este tipo de sistema también es 
conocido con el nombre de levógiro. 


Consideremos un sistema de coordenadas cartesianas ortogonales (x, y, z) con el origen en el punto 
inicial del vector V. El sistema se orienta de tal forma que cumpla la condición de sistema de referencia 
directo. El vector Y queda representado por una terna ordenada de números reales 


v= (vi, v2,v3) 


donde vi, v2 y v3 son respectivamente las proyecciones del segmento rectilíneo correspondiente a 
v, sobre los ejes x,y y z, como se puede observar en la figura . 
Una descripción alternativa del vector Y se obtiene introduciendo los vectores unitarios ej, e2 y e3, 
a lo largo de los ejes x,y y z. Estos vectores unitarios se conocen como los vectores unidad cartesianos. 
También suelen representarse por i, j y K. En función de ellos, el vector Y queda representado por la 
combinación lineal: 
V = vie1 + V2€2 + V3é3 (A.1) 


donde vi, v2 y v3 se conocen como las componentes cartesianas del vector YV. 


lEste producto será definido más adelante en la sección A.4 (Algunos productos de vectores) 
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Figura A.1: Componentes de un vector en un sistema de referencia cartesiano. 


A.4. Algunos productos de vectores 


En esta sección vamos a definir cuatro productos de vectores (escalar, vectorial, mixto y diádi- 
co). Las definiciones atenderán a dos tipos: geométrica, cuando se exprese el producto en función 
del módulo y la dirección del vector; y analítica, cuando se exprese el producto en función de las 
componentes cartesianas del vector. 


A.4.1. Producto escalar o contracción 


El resultado del producto escalar de dos vectores es un escalar (número), de modo que con dicho 
producto se obtiene, en general, un tensor de orden inmediatamente inferior al de los "factores”, de 
ahí que el producto escalar sea conocido también por el nombre de *contracción”. El símbolo que 
representa al producto escalar es un punto centrado (-). 


= Definición Geométrica: 
ú -v= |u||v| cos O (A.2) 


donde 6 es el ángulo que forman los dos vectores. 
Gráficamente, ii - Y representa la proyección de y sobre Y multiplicada por la norma de Y. 


= Definición Analítica. En virtud de que los vectores e1, ez y ea son ortogonales, se cumplen las 
siguientes relaciones: 


€ €; — Óij (A.3) 


donde ő;j se denomina símbolo de kronecker? y se define por: 


1 sii=j 
ó; ¡=> Pia 
0 siisj 
Si expresamos los vectores dí y Y en función de sus componentes cartesianas -ver la ecuación 
(A.1)- y tenemos en cuenta las relaciones anteriores, obtenemos que 


> 


-V= uiv; e) = uiv JÖ; j 


2También llamada matriz unidad de dimensión 3x3 (D o tensor unidad de orden dos. 
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Usando las dos últimas ecuaciones podemos expresar el producto escalar como: 


u-v= uv; (A.4) 


Si tenemos en cuenta que los vectores también pueden escribirse como matrices columna: 


ul 0 0 ui 
i= 0 éi-F uz 82 + 0 83 = u2 (A.5) 
0 0 u3 u3 


el producto escalar de ti y Y será: 


ü- V= uv = (u; uzu3) | va | = uivi (A.6) 


Propiedades del producto escalar 


Algunas propiedades del producto escalar son: 


a (Ai + uV) -w= Mu -w) + u(v -w) siendo A,u € R 


> 


a ¡-Yv=0<= i y Y son ortogonales (si |x| y |V| 4 0). 


A.4.2. Producto vectorial 


El resultado del producto vectorial de dos vectores es otro vector. Dicho producto se simboliza 
con un acento circunflejo centrado (A) o, en ocasiones, también se emplea un aspa (x). El producto 
vectorial, a diferencia del producto escalar, no es un producto de tipo interno porque para su definición 
hace uso de un elemento externo: el símbolo de Levi-Civita*. Sin este elemento no sería posible una 
definición analítica. 


a Definición Geométrica. El producto vectorial de los vectores d y V es un vector de dirección 
perpendicular al plano formado por ambos, de sentido el de avance de un sacacorchos que gira 
del primer vector al segundo por el ángulo menor y cuyo módulo viene dado por: 


la A vI = lüllvi sin 9, (A.7) 
donde 9 es el ángulo que forman los vectores i y Y (ver figura ). 
Gráficamente, [di A V| representa el área de un paralelogramo de lados u y v. 


= Definición Analítica. Para obtener la expresión del producto vectorial en función de las com- 
ponentes debemos proceder de forma análoga a como lo hicimos con el producto escalar. Dado 
que los vectores e, ez y ez son ortogonales y que forman un triedro directo, se cumplen las 
siguientes relaciones: 


3Se definirá más adelante. 
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Figura A.2: Vector resultante de hacer el producto vectorial de los vectores t y V. 


eN éj = Eijkêk (A.8) 


donde la letra griega e representa al llamado símbolo de Levi-Civita*, cuya definición es: 


1 si ijk es una permutación par de 1,2,3. 
€ik=3 —1 siijk es una permutación impar de 1,2,3. 
0 si dos subíndices son iguales. 


Teniendo en cuenta la definición geométrica de producto vectorial y las relaciones anteriores, 
el vector resultante de hacer el producto vectorial de dí sobre Y viene dado según sus componentes 
cartesianas por la ecuación: 

HAVE E¡ ¡ku ¡VLe; (A.9) 

Existe una forma simbólica de realizar el producto vectorial de dos vectores expresados en com- 
ponentes cartesianas. Consiste en efectuar el desarrollo de un determinante un tanto peculiar, que 
como tal no tiene sentido, pero que al ser desarrollado reproduce fielmente las operaciones propias 
del producto vectorial. Sean u = uiel + uz€) + uzé3 y V = v¡€] + Voe? + v3€3, entonces: 


é é; éa 
HxV—| u u uas |= (u2v3 = u3Vo )é1 + (uzvi = u1v3)€2 + (uva = u2v1)é3 (A.10) 
VI Va V3 


con lo que se llega finalmente a la relación (A.9). 


Propiedades del producto vectorial 


El producto vectorial cumple las siguientes propiedades: 
= ¡Aú= —V A ú (propiedad anticonmutativa) 

a ¡AA+ pm) = MUA V) + ul A w) 

a ¡AV=0<= i y Y son paralelos (si [új y |V] 4 0.) 


*También llamado "tensor unitario de orden 3 completamente antisimétrico”. 
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Figura A.3: Representación del paralelepípedo construido sobre los vectores i, Y y w. 


A.4.3. Producto mixto 


El producto mixto se define como el producto escalar de dí por el vector que se obtiene de multi- 
plicar vectorialmente Y y w. El resultado de hacer el producto mixto de tres vectores es un escalar. Se 
puede notar de tres maneras distintas: 


= Definición geométrica. En virtud de las definiciones geométricas del producto escalar y del 
producto vectorial expresadas en las ecuaciones (A.2) y (A.7) respectivamente, obtenemos que: 


> 


t- (VA w) = |ú]|v||w]cos asin B (A.11) 


donde B es el ángulo formado por los vectores Y y W; y 0 es el ángulo que forma el vector ií con 
la normal al plano formado por los vectores Y y w (ver figura ). 


Gráficamente ii - (VA W) representa el volumen de un paralelepípedo de lados u,v y w. 


Definición analítica. Teniendo en cuenta las definiciones de producto escalar y producto vecto- 
rial dadas por las ecuaciones (A.4) y (A.9), podemos expresar el producto mixto como sigue: 


U-(VAW) = (éi) - (Ejkivewiej) = E jujviwi(e; Ej) = E ¡k¡UiVW10;¡ = 
= E€jklU jVkWI 
Como todos los subíndices que aparecen en la última expresión se contraen debido a la suma 


(convenio de Einstein), son simples subíndices mudos, de manera que podemos cambiarles las 
letras. Teniendo en cuenta esto, el producto mixto queda como: 


i- (VA w) — E¡jkUjV jWx: (A.12) 
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De forma análoga a como lo hicimos con el producto vectorial, podemos expresar el producto 
mixto de tres vectores H, Y y W como el desarrollo del determinante 


ú-(vAW) = VI Va V3 (A.13) 
W Ww Wa 


En virtud de esta manera de expresar el producto mixto de tres vectores como el determinante 
formado por las componentes cartesianas de dichos vectores, podemos generalizar y expresar el de- 
terminante de una matriz genérica A en función de los elementos de matriz: 


a a2 413 
lA]=| an 492 az | =E¡jk41¡0203k. (A.14) 
431 432 433 


Pero esta no es más gue una de las posibles formas de expresar el determinante. También se podría 
haber escrito: 
|A| = €; ¡k42543 01x = €ijkd3id1 jak = —Eijkd11d3 jak = tt (A.15) 


Una expresión general gue engloba a todas las posibles combinaciones es: 
|A| = Eimnéijkdtidmjdnk (sin sumar sobre l,m y n.) (A.16) 


Si hacemos la suma sobre l,m y n nos daremos cuenta de que estamos sumando seis veces el 
mismo término, de manera que la ecuación general anterior se puede expresar usando el convenio de 
Einstein de la siguiente forma: 


1 
|A| = g EHk EImn AA ¡man (A.17) 


Propiedades del producto mixto 


El producto mixto tiene las siguientes propiedades, que se deducen a partir de las del producto 
escalar y vectorial. 


a i- (VAW) =1 <= i, Y y w son ortonormales y forman un triedro directo. 


a i- (VAW) = 0 <4 á, Y y w son coplanarios (si lú], |V| y [él 4 0). 


A.4.4. Producto diádico o tensorial 


Se define el producto diádico de los vectores t y Y por su efecto al contraerse con otro vector w a 
través de las relaciones: 


(197) -W=(0-%) (A.18) 


w- (Q7) =(%-1)7 (A.19) 
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Al contraer el producto diádico con un vector w resulta otro vector que tiene la dirección de i si 
se contrae por la derecha -ecuación (A.18)- y de Y si se contrae por la izquierda -ecuación (A.19)-. 


(17) -w= i- w) (A.20) 
w- (197) =(%-1)7 (A.21) 


Teniendo en cuenta que i y Y se pueden escribir en forma matricial, podemos expresar el producto 
diádico como el producto de una matriz columna por una matriz fila: 


u 
9V =uv = | u |(vivzv;)=M (A.22) 
u3 
Siendo la matriz M de la forma M= 


UV] H1V2 UuV3 
M2V1 UVa W2V3 
UzV] U3V2  Uu3V3 


que puede escribirse como 


O = 


0 0 000 
M = uivi 0 0 |+--+u3w3| 000 (A.23) 
0 0 00 1 


© 


Esto puede expresarse como la suma de las componentes de los vectores u y V 
iV = uiv (ei 9 e1) + uive De) +--- +uzv3 (e 9 ez) (A.24) 
A ec; 8 e; se les denomina diadas unitarias. El concepto de diada es una generalización del concepto 
de vector. Las componentes de la diada pueden representarse por los nueves coeficientes de una matriz 
cuadrada 3 x 3. Una generalización de las diadas dará lugar a las poliadas, de 39 componentes. 


Principales propiedades del producto diádico 


a Para a € R se verifica 
(ai) 9 — 19 (av) = a(i 8 Y) (A.25) 


Esta relación es una consecuencia directa de la definición de producto diádico, puesto que las 
direcciones de los vectores no cambian al multiplicarlos por un número real. 


= El producto diádico tiene la propiedad distributiva respecto de la suma, es decir, 
¡SN (V+Ww) =dQV+ 18D W (A.26) 
(+y 9W=d9QW+V9QW (A.27) 
= El producto diádico no es conmutativo: 
US VAVOU (A.28) 


En efecto, 
ud SV= uy j(éi éj) A viu (é Q éj) = VQ U (A.29) 
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Figura A.4: Dos sistemas de referencia cartesianos girados uno respecto al otro pero con un origen 
común. 


A.5. Transformaciones ortogonales de ejes cartesianos 


Una de las transformaciones lineales más sencilla y útil en la interpretación de problemas físicos es 
la ortogonal. Consiste en pasar de un sistema de coordenadas cartesianas a otro sin que haya distorsión 
o cambio de escala. 

Supongamos dos sistemas de referencia cartesianos, como se muestra en la figura donde uno de 
ellos está girado respecto del otro, compartiendo ambos el mismo origen. 

El sistema 1 (sin prima) es el formado por (O,e1,e2,e3), en el que un vector cualquiera se expresa 
como combinación lineal de los vectores de esta base como: 


V= vie; (A.30) 


4 4 a >>>! : 
mientras que el sistema 2 (con prima) es el formado por ([O.e1 ,e2 ,e3 ), en el que un vector cualquiera 
se expresa como combinación lineal de estos vectores como: 


ej (A.31) 


Para poder determinar las coordenadas del sistema 1 a partir de las coordenadas del sistema 2, 
necesitamos conocer las componentes de los vectores básicos de un sistema en función de la base del 
y >! $ é š $ 
otro. Por ejemplo, el vector e; de la base 2 (con prima) referido a la base 1 (sin prima) se expresa 
como: 


= z$ 
ei = Aij; (A.32) 
donde A;j son los cosenos directores. Usando la ecuación (A.30) podemos expresar un vector genérico 
V como: 
V= vé; 
pero haciendo uso de la ecuación (A.31) podemos expresar el mismo vector haciendo referencia 
a otra base: 
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y como se trata del mismo vector, v, podemos igualar las dos expresiones anteriores: 


Vie; = vie; (A.33) 


Si ahora hacemos uso de la ecuación (A.32) podemos escribir 
v¡é = VA jkêk 
Multiplicando escalarmente por e; a ambos lados, se obtiene 
Vi = viji (A.34) 


gue nos da la relación de las coordenadas de un vector genérico en un sistema con las correspondientes 
a otro sistema girado respecto a éste. 
La transformación inversa se otiene multiplicando escalarmente por e” y la expresión A.33, siendo 
el resultado 
vi = A jivi (A.35) 
Como ambas bases son ortonormales, se cumple gue: 
ejej = Óij (A.36) 
El Ej = Óij (A.37) 
Sustituyendo la ecuación (A.32) en la ecuación (A.37) obtenemos 


Si — ej € = (Mixer): (Are) = Mukk jiëk e) = MIRA jÒ = MA jk 


de donde deducimos gue 
dir A jk = Oj (A.38) 


Esas seis ecuaciones ponen de manifiesto que de los nueve cosenos directores (A);;, sólo tres son 
independientes. Esta ecuación puede expresarse también como 


dir (Ag) = Sij (A.39) 


que nos recuerda a la expresión en componentes del producto de dos matrices. El resultado, según 
corresponde, será una matriz de orden 3 x 3. La delta de Kronecker ó;; puede entenderse como la 
expresión en componentes de la matriz unidad de orden 3 x 3, en la que los elementos de la diagonal 
(i = j) valen 1 y el resto (i sé j) vale 0. Por tanto, la ecuación (A.39) se escribe en forma matricial 
como: 


AA =I (A.40) 
Observando esta última ecuación bajo la perspectiva de las propiedades matriciales, se hace evidente 
que la matriz A! es la inversa de A, por tanto se cumple que: 


A =A! (A.41) 


En componentes 
Ap =Aji (A.42) 


1 


En definitiva, podemos concluir que la matriz A es ortogonal. 
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Por dicha ortogonalidad , podemos escribir la transformada inversa: 
El = Ne (A.43) 
E =A E = Aji)! (A.44) 


Según la transformación sufrida por una magnitud tras una rotación, podemos introducir las sigu- 
ientes definiciones: Una magnitud V cuyas componentes, tras una rotación del sistema de coorde- 
nadas, se transforman de la forma 

1 
v= Aijvj (A.45) 


es un tensor de rango uno. 
Una magnitud T cuyas componentes, tras una rotación del sistema de coordenadas, se transforman 
de la forma 
1 
T;; = My Ta (A.46) 


1 


es un tensor de rango dos. El tensor T puede expresarse en función de las diadas unitarias 
T = Tijés DÉ; (A.47) 


Por generalización, un tensor T de rango n es una magnitud que tiene 3” componentes y que bajo 
una transformación ortogonal sufre el siguiente cambio: 


Ti. an = Maji Majo i MINJNT 89... jN (A.48) 
T podemos escribirlo en función de las poliadas: 
T = Tig. Néit 9 ép a : Q Cin (A.49) 


Las diadas y poliadas unitarias constituyen la base de los tensores de rango dos y N respectivamente. 


A.6. Invariantes 


La mayor parte de las condiciones de simetría espacial recaen sobre las magnitudes gue aparecen 
en las leyes físicas, pues deben poder expresar toda la variedad de sistemas y de leyes gue existen. 

En esta sección nos centraremos en la invariancia frente a cambios de sistemas de referencia 
ortogonales de los tensores y de los productos escalar, vectorial y diádico. 

Consideremos dos sistemas de referencia con el mismo origen pero girado uno respecto al otro. 
Vamos a realizar una transformación ortogonal para poner de manifiesto gue los productos escalar, 
vectorial y diádico son invariantes frente a cambios de sistemas de referencia. Vamos a denotar por u; 
y v; alas componentes de los vectores dí y Ven un sistema de referencia, y por u, y v; a las componentes 
en el otro sistema de referencia. 


A.6.1. Invariancia del producto escalar 


Vamos a demostrar que el producto escalar es invariante bajo cambios de coordenadas ortogonales. 
Sean i y V dos vectores en las dos bases dadas a continuación 


ej (A.50) 


v= Vie; —vjej (A.51) 
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el producto escalar de ti por Y será: 


U-V= Uv; = (A jruk) (A jivi) = A jÀ ¡ru! (A.52) 
donde se han usado las relaciones 
| = Ajik (A.53) 
= hj (A.54) 
dado gue 
A jkA ji = Ôr (A.55) 


Podemos escribir la expresión (A.52) de la forma: 
U-V=u vi = Ök UkVI = UKVk (A.56) 


por lo gue gueda demostrada la independencia con respecto al sistema de referencia. 


A.6.2. Invariancia del producto vectorial 


El producto vectorial de los vectores t y V en las distintas bases se expresa a continuación: 


Ü NV = Ej ¡gu jVxéi (A.57) 
úNv = Ej ¡gl vE (A.58) 
teniendo en cuenta las relaciones (A.57) y (A.58) junto con las expresiones 
€! = Mimém (A.59) 
Eijk = Pit jmAknEmn (A.60) 
podemos escribir 
úNv= Es ULVE e = (AA jmMin€lmn ) (A nun) (Arvi) (A imém) (A.61) 


y si recordamos la definición de la delta de Kronecker la ecuación anterior quedará de la forma 
Sr> / ikai > > 
UNV = E; ¡pu Ve = ÒimÒmh Òn EImnUhV Em = ElmnUmVnël (A.62) 


Oueda así demostrada la invariancia del producto vectorial. 
La transformación de €;/x en las transformaciones ortogonales se comentará en la sección de ten- 
sores isótropos. 


A.6.3. Invariancia del producto diádico 
El producto diádico en dos sistemas de referencia ortogonales vendrá dado por las expresiones 
úgv = ujvj (e, éj) (A.63) 
ügyv — uv (e 8 E) (A.64) 
por tanto, si partimos de la expresión (A.64) y recordando las expresiones (A.53), (A.54) y (A.56) 
obtenemos la expresión: 


úgv — uv (e D éj) = (Aitte) (A jivi) (Nimm) Q (Ajnén) (A.65) 
si usamos ahora la relación (A.55) en la expresión anterior quedará demostrada la invariancia del 


producto diádico: 


e És >lo>! = S = = 
US V — ug (ë; DE ) = ÖrmörmukvIi ém Q én = ULVIE O éj (A.66) 


i 
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A.6.4. Tensores isótropos 


Un tensor es isótropo cuando sus componentes no varían bajo una rotación, es decir, es invariante 
frente a rotaciones. 


m Tensores de orden cero. Los tensores isótropos de orden cero son los escalares. 
= Tensores de orden uno. No hay tensores isótropos de orden uno. 
= Tensores de orden dos. Los tensores isótropos de orden dos son de la forma A;; = Qð; j, Siendo 


Oo. una constante. 


La transformación de 9; ¡ en una rotación de ejes será: 


1 


87 a Mik A 51611 = Ag) jk = Ő; j (A.67) 


= Tensores de orden tres. Los tensores de orden tres son de la forma Ay; = Be; jk, donde B es una 
constante. 


Un ejemplo de tensor isótropo de orden tres es el símbolo de Levi-Civita, £;jx, que como ya 
hemos visto es usado en la definición del producto vectorial. Basándonos en la invariancia del 
producto vectorial, podemos concluir afirmando que dicho símbolo tiene que ser forzosamente 
invariante. 


Veamos ahora una demostración más detallada. Queremos evidenciar que la forma de €;jx es la 
misma que la de AA ¡mArn€lmn- 


La expresión que relaciona al tensor de Levi-Civita en dos sistemas de referencia distintos es: 
Eijk SS A (A.68) 

siendo A/A jmMtn€lmn el desarrollo del determinante de la matriz A. 

Según habíamos definido anteriormente, la matriz A cumple las propiedades siguientes: 

0 siis 7. 

0 sii=j=k. 

1 para permutaciones pares de 1,2,3. 

1 para permutaciones impares de 1,2,3. 


/ P 
Por tanto, e; jk Será 


0 si hay subíndices repetidos. 
E; jk = 1 para permutaciones pares de 1,2,3. 
—1 para permutaciones impares de 1,2,3. 


que es precisamente la forma que tiene €;jx. Por tanto, vemos que £; ix tiene la misma forma que 
Eijk, lo que demuestra la invariancia. 


a Tensores de orden cuatro. Los tensores de orden cuatro son de la forma: 


Ai jkl = ASi jÒ + VO ¡x0 ¡1 + Y9 ¡10 jk (A.69) 
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A.7. Producto interior de tensores (contracción) 


Sean T y R dos tensores de orden dos 
T=) 
R= Rijei 8 €; 
se define el producto interno o contracción de T y R como 
T-R = (T; 8 E (Ruk D €1) = Tjele ex) O Ruei = 
= TijRaei 9 ejőkj = TR je; O €) (A.70) 


Con el producto interno lo gue conseguimos es disminuir el orden total en dos unidades. 
También podemos considerar la doble contracción de T con R 


T:R=T5R; (A.71) 


cuyo resultado es un escalar. 
Si contraemos un tensor de segundo orden consigo mismo obtenemos 


T:T=T; (4.72) 


que es la traza del tensor 7. Esta cantidad, como ya hemos comentado es un escalar verdadero y por 
tanto un invariante. 

Veamos la actuación de un tensor de segundo orden sobre un vector: 

Sea v el vector resultante de la contracción del tensor T sobre el vector i. 


V =T-ú=T;¡0, 9 €; urek — T;jux(e; 9 €) ex = 
= Tijuxei(e; ex) = Tijure;Ó jk = Tiju jê; = Vie; (A.73) 
La contracción también podríamos haberla hecho por la izquierda: 
v =ú-T = uxrT;jek(e; 9 Ej) = uxTi(ex " E)E = 
= ur Tökre j = uilijej = Vjej (A.74) 


A.8. Producto exterior de tensores 


Sean T y R dos tensores de orden dos 
T =T¡¡6, 9 €; 
R=R;j8¡ 3 €; 
Definimos S como el producto tensorial de T y R de este modo 
S=T8R=T;¡¡Rue D Ej De 8 € = Siju Ei De VEB ei (A.75) 


El rango de S viene dado por la suma de los órdenes de R y T, por tanto, en este caso S es un tensor 
de orden cuatro. 
La generalización a orden N es la siguiente: 


T = Tig. iN€1 2 €2 2 : DEN (A.76) 


R=R ¡N+1jN+2...¡2NÉN+1 Q Én+2 Q` O EN (A.77) 
S=T QR = Ta. iwRjN+1..¡2né1 69: 60 én Vén+1 69: S ÉN (A.78) 
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A.9. Autoanálisis de tensores de segundo orden 
Sean ú y T un vector y un tensor de segundo orden respectivamente, 
T =T,¡é¡ 9 €; 
Ú = ukêk (A.79) 


Se denomina autoanálisis del tensor T al estudio de las soluciones de la ecuación 
T -ü= An (A.80) 


donde ú se denomina autovector y À autovalor. Es decir, es la búsqueda de direcciones del espacio 
en las que T actúa simplemente como un escalar. Para resolver el problema planteado es conveniente 
reescribir la ecuación anterior pasando todos los términos al primer miembro de la ecuación: 


(T-AD)-ú=0 (A.81) 


La condición necesaria y suficiente para que este sistema homogéneo tenga solución distinta de la 
trivial (zi = 0) es que el determinante de sus coeficientes sea cero: 


det(T —M) =0 (4.82) 


Si desarrollamos este determinante obtenemos un polinomio cúbico denominado polinomio carac- 
terístico. La ecuación característica tiene la forma 


E=BLA+LM-M=0 (A.83) 


De esta ecuación extraeremos tres soluciones, A"), 12) y 10) que podrán ser, en general, una real y las 
otras dos complejas conjugadas. A partir de estos autovalores podremos calcular las direcciones prin- 
cipales o autovectores ti. Esto se entenderá mejor viendo el ejemplo del final de la sección Invariantes 
del tensor de segundo orden. 

En el caso real y simétrico los tres autovalores serán reales y los autovectores correspondientes a 
autovalores distintos serán ortogonales. 

Los autovectores los podemos elegir de forma que constituyan una base ortonormal, ív 0) 70) 46) ). 
Si se cumple 

y. (10) AO) =1 


entonces (50,50) 78), forman una base ortonormal directa y son el sistema de ejes principales del 


tensor T. En el sistema de ejes principales el tensor es diagonal. Cuando nos refiramos al tensor T en 
el sistema de ejes principales lo denotaremos por T*. 
Lo que hemos conseguido es una nueva base donde el tensor es diagonal. Vamos a demostrarlo: 
Si se tiene 
T —T;jé 9 é; (A.84) 


entonces la componente 7; ; puede obtenerse como 
Tij = ĉi T "éj = éj" (Taek 9 é1) "éj = Tu (é ék) (ér: éj) = Tu101 01; = Tij (A.85) 
Las componentes de T en la base de los ejes principales, e; " ortonormal, son: 


T¡*=8*-T éjt (A.86) 
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Pero como éj" es el autovector correspondiente al autovalor 20), entonces: 
Tt = gae * = Agy. éjt = DE; (A.87) 
En virtud de los últimos resultados, podemos escribir 
T* = APS g 808 *=1V8 80841202412 eg (A.88) 


en los ejes principales. 
Hemos demostrado de esta manera, gue el tensor en los ejes principales lo podemos escribir como: 


Tá = Mk Ta = A’ ô; (A.89) 
AD 0 0 
T*=| 0 A 0 (A.90) 
0 0 28 


A.10. Invariantes del tensor de segundo orden 


Sea T un tensor de segundo orden, T = 7; je; 2 ej. En el espacio tridimensional, de todo tensor de 
segundo orden se deducen tres invariantes, que se conocen como invariantes canónicos: 


a Invariante lineal. Este invariante coincide con el determinante de T. El determinante es un 
invariante, como ya hemos visto, ya que es un escalar verdadero. Por tanto el valor del determi- 
nante de T es el mismo en todos los sistemas de referencia. 


1 
h = det(T) = ¿eipemTaTmTin = 1101018 (A.91) 


= Invariante cuadrático. Mediante un argumento análogo al anterior vemos que esta cantidad es 
un invariante. 


b -| hí Tp la 113 Es 123 |_ 
Tu 122 Bi 133 Ta; 133 
1 
= z TiTi —T;jTj) = 
= ADA +A 101.6) (A.92) 


a Invariante cúbico. Esta cantidad es un invariante ya que coincide con el valor de la traza de T, 
que es un escalar verdadero. 


B=tr(T) = Ti = AD 4 AC) EAG) (A.93) 


Ejemplo: Sea el tensor A;j, cuyas componentes son los elementos de la matriz 


1 —4 6 
Aj=| -4 -3 -2 (A.94) 
ő. 2 2 
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Oueremos reducirlo a la forma diagonal por una transformación ortogonal. La ecuación característica 
o secular es 
l-o -4 6 


4 3-06 2 |=-0%+630+162=0 (A.95) 
6 —2 2—o 
cuyas raices son 
01 = —3 
02 = 9 
03 = —6 


Por tanto, la forma diagonal buscada es 
—3 0 0 
Aij — 0 9 0 (A.96) 
0 0 


La transformación ortogonal gue nos permite pasar de una base a otra se obtiene resolviendo los 
sistemas: 


=a Para ©] = —3: 
40 —4B+6y=0 
—-4a—27—0 
6a— 28 -£57- 0 
cuya solución en función de al es 
y= —204 
B = —2a 


De donde se obtiene el autovector (ya normalizado): 


pl 1 
ui = 3 (—1 252) 
a Para 62 = 9: 
—8a — 4B +6y=0 
—4a — 128 -2y=0 
6a —2B —7Y— 0 
cuya solución en función de B es 
a — —2B 
Y— —2B 
De donde se obtiene el autovector (ya normalizado): 
M 1 
U = 321,2) 
a Para 63 = —6: 
70. —4B+6y=0 
—40+3B-2y=0 
60 —2B+8y=0 


174 APÉNDICE A. NOCIONES BÁSICAS SOBRE CÁLCULO TENSORIAL 


cuya solución en función de yes 


Chris: 
y=% 
De donde se obtiene el autovector (ya normalizado): 
1 1 
lia = ——- (—13,5,11) 


v315 


A.11. Tensores en coordenadas curvilíneas 


Podemos definir los tensores de modo más general basándonos en una transformación arbitraria. 
Para ello consideramos dos conjuntos de variables independientes, x; y xi. Entre estos dos conjun- 
tos de variables definimos la siguiente transformación: 
/, /, 


x; =X;(x;) (A.97) 


Ahora bien, para que esta transformación esté definida correctamente es necesario que sea invert- 
ible. Es decir, debe ser posible escribir tanto la relación (A.97) como: 


xi =xi(x;) (A.98) 


Para que esto ocurra es necesario y suficiente que el jacobiano de la transformación sea diferente 
de cero en todos los puntos.Se llama jacobiano al determinante de la matriz jacobiana. Por su parte, 
la matriz jacobiana no era más que la matriz que relacionaba unas variables con otras a través de 
las derivadas parciales. Por tanto, la condición necesaria y suficiente para que la transformación sea 
invertible es: 


pa peca EST] (A.99) 


De forma más explícita podemos escribir la expresión (A.99) como: 


o... d 
Ox, dx, ox, 
Ox dx Ox 
fe 2 "(120 
dn Mn... An 
Ox, dx, ox, 


Visto esto, cabe cuestionarse la manera de expresar de forma más genérica los escalares, los vec- 
tores, los tensores, etc. 

En primer lugar construyamos la magnitud que vamos a usar como referencia para definir un 
vector: 


1 
, Ox ; 
dx; = dx; (A.100) 
Ox J 
Cualesquiera tres magnitudes que en la transformación dada por (A.97) se transforman como la 
cantidad dx; se dice que forman un vector contracovariante. Pues bien, a las componentes de este tipo 
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de vectores se les colocan superíndices. Por tanto, debemos introducir una nueva notación. Así pues, 
la transformación general dada por (A.100) la escribimos ahora como: 


HA 


xi = x (a) (A.101) 


Por consiguiente, la transformación entre los vectores contravariantes se define como: 


7 


a a, 
A! = — AJ (A.102) 
Ox! 
No obstante, como es evidente, no todas las magnitudes de tres componentes se transforman de esta 
manera. Por ejemplo, el gradiente de un escalar presenta una transformación distinta. Veámoslo: 
Suele definirse el gradiente de un escalar,A, como: 


JA 
B; = — (A.103) 
dxi 
La ley de transformación en un cambio de coordenadas será: 
JA ÐA Oxi 
E = — = — : — A.1 4 
2 axi Oxi axi Coinn 
y, si cambiamos el orden, nos gueda: 
; JA Ox] JA ox 
B; siet E (A.105) 


ax ax axi Oxi 7 


Pues bien, cualquier magnitud A que se transforma según la relación (A.105) es un vector covari- 
ante. En este caso, para nombrar las componentes de dicho vector se emplean subíndices. 

Ahora bien, podemos pensar que los vectores contravariantes y covariantes están relacionados ya 
que cualquier vector puede ser expresado tanto en componentes contravariantes como en componentes 
covariantes. El tensor que relaciona ambos "tipos"de componentes es el llamado tensor métrico. 

El tensor métrico no es más que el conjunto de funciones que definen la métrica del espacio, es 
decir, la distancia entre dos puntos del espacio. En el caso del espacio tridimensional se tendrán seis 
funciones independientes.Podemos hacernos las preguntas ¿Cómo definimos la métrica de un espacio 
cualquiera? ¿Cómo medimos distancias en cualquier sistema de referencia? Si tenemos dos puntos 
infinitamente próximos, la distancia entre ellos viene dada por la ecuación: 


ds? = gijdx dx! (A.106) 
Y esta expresión recibe el nombre de primera forma fundamental, donde: 
Sij = gij(x") (A.107) 


es el tensor métrico, el cual es simétrico. 


Ejemplo: Veamos el tensor métrico en el caso de las coordenadas cartesianas en el plano. 
ds? = dx" +dy? 


donde: 
81 =1 
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gn =1 
gn =0 
y este tipo de métrica se puede resumir en: 
gij = Óij (A.108) 


Todas las métricas cuyo tensor métrico se hace cero cuando los subíndices son distintos entre sí 
reciben el nombre de métricas ortogonales. 


Ejemplo: Veamos las coordenadas polares en el plano. 


ds? = dr’ + r"de? (A.109) 
donde 
8r=1 
800 = r 
gro=0 


La distancia entre dos puntos no viene dada por el teorema de Pitágoras. 


Volvamos al concepto de tensor métrico. Siempre debemos tener presente que el tensor métrico 
nos informa de cómo se mide en el espacio en el que se trabaje. Por otra parte recordemos que hemos 
introducido el tensor métrico para encontrar la relación existente entre los vectores covariantes y 
contravariantes. Veamos, pues, tal relación: 


A¡=gij4! (A.110) 


Observando la expresión anterior se ve que cuando se contrae un tensor contravariante con un tensor 
métrico, se obtiene un tensor covariante. Nótese que el tensor métrico tiene subíndices en lugar de 
superíndices. Esto es debido a que es covariante. Volvamos a la primera forma fundamental. Según lo 
que acabamos de ver: 

ds? = gy¡dx dx! = dx (gijdx!) = dx dx; (A.111) 


Según esta última ecuación, la distancia entre dos puntos no varía; es decir, es un escalar verdadero. 
En general tenemos que: 


a Si contraemos un vector covariante con un vector contracovariante, obtenemos un escalar ver- 
dadero. 


= Si contraemos un vector covariante con otro vector covariante, no se obtiene un escalar. 


= Si contraemos un vector contracovariante con otro vector contracovariante, no se tiene un es- 
calar. 
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Ahora bien, según todo lo gue hemos visto hasta ahora, la operación inversa de (A.110) resulta ser: 
AJ =gUA, (A.112) 


Es evidente que g' son las componentes contravariantes del tensor métrico (de ahí que se usen su- 
períndices en lugar de subíndices).Por definición: 


g” = (gi) (A.113) 


Consecuentemente: 
818 = aí (A.114) 


Como vemos, la expresión anterior no es más que una nueva definición de la delta de Krönecker. 
¿Qué ocurre si estamos en coordenadas cartesianas? Como las "g"son la "6", las componentes covari- 
antes coinciden con las componentes contravariantes (es decir, son las mismas). Esto es lo que hace 
que empleemos siempre subíndices para nombrar a las componentes de los vectores sin preocuparnos 
de si son contravariantes o covariantes. 

Llegados a este punto debemos preguntarnos por las transformaciones entre tensores. ¿Cómo 
definimos un tensor de segundo orden dos veces covariante? ¿Cómo lo definimos si es dos veces 
contravariante? ¿Cómo lo definimos si es mixto: una vez contravariante y otra covariante? De acuerdo 
a la transformación de coordenadas que hemos visto, se definen como sigue: 


= Sies dos veces contravariante: 


1. xix! 
yu = ETA A.115 
oxk dx! ( ) 
u Sies dos veces covariante: . 
/ Ox dx! 
ii ( ) 


a Si es mixto: covariante y contravariante: 


AS 
ij Ox ox! l 


real (A.117) 


Ahora bien, ¿cómo pasamos de contravariante a covariante, es decir, cómo se suben y se bajan 
índices? Con el tensor métrico. Por tanto: 
_ kl 
Tij = 8x8 jT (A.118) 
T} = gT” (A.119) 


Esto contiene el caso de las transformaciones de las coordenadas ortogonales: 


g= N (A.120) 
y de aquí se saca que: 
dx, 
Aij = — A.121 
o ( ) 


y de 
xj =2x; (A.122) 
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Ox; 
TE z (A.123) 
i 
entonces a i a 
X; Xj 
Mij = T = Sr (A.124) 
J i 
y, por tanto 
Ty; = Mala =T” (A.125) 


de aquí que no haya que hacer diferenciación entre componentes covariantes y contravariantes en un 
sistema cartesiano: j Se transforman de la misma manera! Así pues, de aquí se deduce también el que 
no se usen unas veces superíndices y otras subíndices. No obstante, hay materias de la física en la 
que esto no se puede obviar. Es el caso de la relatividad, en la que sí es importante diferenciar entre 
vectores covariantes y contravariantes, siempre que las coordenadas sean curvilíneas. 


A.12. Diferenciación de tensores cartesianos 


Supongamos ahora que tenemos un campo tensorial representado por un tensor de segundo orden 
siendo las componentes del tensor funciones de las coordenadas. Así pues, nuestro tensor resulta ser: 
A¡¡(xx). Pues bien, vamos a construir la siguiente magnitud: 


JA; j 
S (A.126) 
OXk 

Vamos a ver la actuación del operador nabla, v, en coordenadas cartesianas sobre un tensor de se- 
gundo orden. La actuación de nabla sobre un tensor sin contraer se denomina gradiente. La actuación 
de nabla sobre un tensor contrayendo se denomina divergencia. 

Para referirnos al gradiente en varias dimensiones se usa el operador nabla: V. Por su parte, el 
gradiente se puede representar también por Ve o por V. El que detrás de nabla no aparezca nada o 
aparezca el símbolo ® significa que "no se contrae". 

Veamos cómo se transforma la siguiente ecuación en un cambio de coordenadas ortogonales: 


JA; 
Bij = 2 A.127 
jk — P ( ) 
Podemos escribir 
JA; ; 9 MA ¡mA m OXn OA im 
ijk = A = A e is = = Nil jmen = m 

X, X dx, Xn 

Ail À men Bimn (A.128) 


y esto no es más gue el gradiente. 

Si observamos la expresión anterior, vemos gue el gradiente ha hecho gue el tensor nos aumente 
en un grado; es decir, teníamos un tensor de segundo orden y al hacer el gradiente hemos obtenido un 
tensor de tercer orden. Así pues, en general, al hacer el gradiente de un tensor de orden N se obtiene 
un tensor de orden N +1. 

El resultado anterior no es correcto en el caso de tensores no cartesianos. 
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Antes de continuar, consideremos la siguiente notación: 


9Aij =Ajok (A.129) 
OXk 
como vemos, cuando tengamos uno o más subíndices separados del resto de subíndices por una coma, 
tendremos que derivar la expresión respecto a la componente que tenga esos subíndices. En definitiva, 
una coma dentro de un subíndice implica una derivada. 
Volvamos al operador nabla. El operador nabla’ se define como: 


y_% (A.130) 
oxi 


¿Cuáles son las operaciones que se pueden realizar con este operador sobre cualquier magnitud? 
Estas operaciones dependen del tipo de magnitud sobre la que actúen. Así pues, tenemos: 


= El operador nabla actuando sobre un escalar. 


En este caso el operador no puede actuar más que de una manera: como gradiente. Por tanto, si 
tenemos un escalar A, su gradiente vendrá dado por VA. Dicho gradiente es un vector: 


Vie (A.131) 
oxi 


= El operador nabla actuando sobre un vector. 


En este caso el operador nabla puede actuar sin contraer o contrayendo. 


Por tanto, si tenemos un vector B = B; , nos podemos encontrar con: 


f >_>- 0B 
Gradiente de B = VB = — (A.132) 
Ox j 
donde vemos gue no se contrae sobre ningún índice puesto gue no hay índices repetidos. Consecuente- 
mente, el gradiente de un vector es un tensor de segundo orden. 


> > > OB; 
Divergencia de B = V-B = F (A.133) 
Xi 
La divergencia de un vector es un escalar. 
e a 9B 
Rotacional d B = V A B = eina (A.134) 
Xj 


donde vemos que para obtener un vector hemos tenido que introducir un tensor de tercer orden. 


Ejemplo: Demostración en coordenadas cartesianas de que se verifica la siguiente identidad: 


VA(fA) = (FAAA +IVA) (A.135) 


5Este operador siempre en un vector, se le ponga o no la "flechita"encima. 
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Escribamos el rotacional de fA en componentes teniendo en cuenta gue B-VA ( fA), por tanto: 
B; = E: (SAL), j= Ei (FAL ¡ES ¡Ax) = fEijkAk j +Eijk f ¡Ar (A.136) 
Esta ecuación la podemos expresar vectorialmente: 
VAA = f(VAA) + (VA) AA (A.137) 


Si comparamos (A.137) con (A.135) vemos gue son idénticas salvo en el orden de los sumandos del 
segundo miembro. Como el orden de la suma (en este caso) no afecta al resultado, gueda demostrada 
la identidad (A.135) gue es lo gue se nos pedía. 


Ejemplo: Prueba en coordenadas cartesianas de gue es cierta la siguiente identidad: 
(A-V)f=A Vf (A.138) 
Expresemos el primer miembro de la igualdad en componentes: 


asr (a2) = CAE (A.139) 


Ejemplo: Demostración en coordenadas cartesianas de que es cierta la siguiente identidad: 


VMANB)=A(V-B)-B(V-A)+(B-V)A—(A-V)B (A.140) 


Al igual que en los casos anteriores vamos a trabajar con componentes pero antes multiplicamos 
el rotacional del producto vectorial de los vectores A y B por Ép: 


Ép: [VA (4 A B)] =8parlErjiA jBi] = EpgkEkjiA ¡Birg +EpgkEkjiA ¡ra Bi (A.141) 
Por otra parte, recordando la definición de la delta de Krönecker, podemos escribir: 


Ep: [VA (4^ B)] = (Šp jõgi — Őpiógj JA jBisq 1 (Őpjögi — Öpiðq;)A jq B 
SA BA Beat Apg B Apah (A.142) 


Y reescribiendo la expresión anterior en forma vectorial tenemos: 


> > > 


VAMANB)=A(V-B)-(V-A)B+(B.V)A—B(A-V) (A.143) 


Observando la ecuación anterior vemos que es igual a la ecuación (A.140) salvo en el orden de los 
términos. Dicho orden, en este caso, no afecta el resultado. Así pues, queda demostrada la expresión 
(A.140). 


Ejemplo: Prueba en coordenadas cartesianas de la veracidad de la expresión: 


VA(VAA) =V(V-A) -VA (A.144) 
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Vamos a seguir los mismos pasos gue en el ejemplo anterior: 


Ép: (V AA) = £ijkÅk,j (A.145) 
por lo tanto: 
Ép: [V A (Ý AA)] = EpqrlErjkA j, j lq = EpqrErjkAk pk A pogg 
= (Sp; Sar SS ÖpkÖgj Ar ją Fa 
= Aropk —Aprgg (A.146) 


y escribiendo la ecuación (A.146) en forma vectorial queda: 


>> 


VA(VAA) =V(V-A) — V?A (A.147) 


que es la expresión que queríamos demostrar. 


A.13. Teorema de Hamilton-Cayley 


Toda matriz cuadrada satisface su ecuación característica. Su demostración es complicada y 
tediosa, por tanto no la vamos a hacer. Sin embargo, vamos a ver qué significa. Para ello vamos a ver 
un caso particular. 

Supongamos que tenemos una matriz cuadrada cualquiera de orden N. Llamemos a dicha matriz 
A. Su polinomio característico es de orden N y viene dado por: 


det(A — xI) = |A — xl| (A.148) 
Por tanto, su ecuación característica es: 
f(x) =|4—x1| =0 (A.149) 


donde A es la matriz cuadrada dada, / es la matriz identidad y x es un número que permite que se 
satisfaga la ecuación. 

El teorema de Hamilton-Cayley dice que si sustituimos el número x por la matriz cuadrada A, se 
sigue satisfaciendo la ecuación característica, esto es: 


F(A) =0y (A.150) 
donde Oy representa a la matriz cuadrada nula de orden N. Expresado en forma matricial: 
0 0 0 
0 0 0 
0 0 0 


bI— AS IA? — A? = (A.151) 


Ejemplo: Supongamos gue tenemos la siguiente matri7 cuadrada 2 x 2: 


1 3 
a-( , J (A.152) 
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En primer lugar, calculemos su ecuación característica: 


IA —x1| = 


jé A SE ET 0) (A.153) 


4 5-x 


La ecuación característica será por tanto: 
x—6x—7—0 (A.154) 


Es evidente gue hay varios valores de x gue satisfacen esta ecuación. Dichos valores se obtienen sin 
más gue resolver la ecuación. Ahora bien, el teorema de Hamilton-Cayley dice gue la matriz cuadrada 
A satisface dicha ecuación, es decir, gue se cumple: 


TAON MÉ 
(0) 1): 


1 
4 
13 18 70 
= (3 7 )- (a 30 )- ( 3 )=0 st 


Así hemos comprobado el teorema para nuestro caso particular. 


Veamos si es cierto: 


A.13.1. Corolario del teorema de Hamilton-Cayley 


Cualquier potencia de una matriz cuadrada A (m x m) se puede expresar como una combinación 
lineal de las matrices I, A, ..., A"! siendo A” la potencia n-ésima de A, donde los coeficientes de la 
combinación son función de los invariantes de A. 

Para aclarar ideas, veamos un caso particular: 

Supongamos que A es una matriz 3 x 3 y que buscamos A. 

Como A satisface la ecuación característica, ecuación (A.151), podemos de ella despejar A? como 


A? = hI -hA + hA’ (A.157) 


Efectivamente, hemos podido expresar la potencia A? en función de I, A y A?. Pero, ¿qué ocurre si 
nos piden A4? ¿Podremos expresarla en función de /, A y A? o necesitaremos también expresarla en 
función de A"? Según el corolario no será necesario A?, ya que 


At =LA-LA? +14 =bA- hA’ +h (b1—LA+1,4?) = 
(1B)1+(B-Kb)A+ (Ñ —Lb) A? (A.158) 


Como podemos ver, efectivamente la máxima potencia es m — 1; en nuestro caso 3 — 1 = 2. De este 
caso particular podemos deducir la expresión general: 


= bml +dm1A + bm 4 +--+ b14”! (A.159) 


donde los b; son funciones de los invariantes de A. 
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A.13.2. Función de una matriz 


Para definir la función de una matriz es necesario conocer el desarrollo en serie de potencias 
de Taylor de esa función para una variable "habitual". En dicho desarrollo se sustituye la variable 
"habitual"por la matriz. Para verlo más claro, hagamos un ejemplo concreto: 

Supongamos gue tenemos el seno de una matriz A cualguiera de dimensión m x m, es decir: 


f(A) =sin(A) (A.160) 


La pregunta que se nos viene a la cabeza es ¿qué significado tiene esta expresión? Pues bien, lo que 
tenemos que hacer es remitirnos al desarrollo en serie de potencias de Taylor del seno. Es decir: 


pr sl 


MASA O 3? 
Una vez conocido el desarrollo oportuno, para definir el seno de la matriz basta con sustituir en el 
desarrollo de Taylor la x por la matriz correspondiente: 


Tar (A.161) 


AZ AS A 

MASA 

Una propiedad interesante es que, a pesar de que el desarrollo en serie de potencias de Taylor del 

seno de x tiene infinitos términos, el desarrollo del seno de la matriz A no tiene infinitos términos 

puesto que, en virtud del corolario del teorema de Hamilton-Cayley, podemos expresar cualquier 

potencia de A como combinación lineal de matrices, siendo A”” ! la matriz de mayor potencia si la 
matriz A es de dimensión m x m. Para el caso m = 3, se tiene: 


hm (A.162) 


sinA = caí + c2A +c¡A? (A.163) 


Los coeficientes que acompañan a las potencias de la matriz en la combinación lineal, aunque son 
un número finito, pueden estar formados por una suma infinita de términos; es decir c1, co y ca (en 
nuestro caso) pueden ser una suma de infinitos términos, funciones de los invariantes de A. 
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En esta sección vamos a estudiar dos formas distintas de descomponer un tensor de segundo orden 
en dos sumandos con unas características determinadas. Por un lado lo descompondremos de manera 
tal que uno de los sumandos tenga como peculiaridad su isotropía. 

Mediante la segunda forma de descomposición obtendremos un sumando perfectamente simétrico 
y otro antisimétrico. 


A.14.1. Partes esférica y desviador 


Se denomina tensor esférico a un tensor isótropo, que puede descomponerse en dos partes,A = 


Az +Ap, donde la parte esférica hace referencia a la componente isótropa, es decir: 
1 
Ag= 3 Tr(A)I (A.164) 


y la parte desviatoria o desviador se corresponde con la parte restante: 


1 
Ap-A-—s Tr(A)I (A.165) 
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Az tiene una componente independiente, mientras gue Ap tiene, por tanto, ocho componentes 
independientes. La relación gue cumplen las componentes del desviador es gue su traza es nula, es 


decir: 
Tr(Ap) =0 (A.166) 


Esta propiedad se demuestra de forma trivial a partir de las propias definiciones: 


Tr(Ap)=Tr(4) —Tr(A) Tr(1) = Tr(A) —Tr(4)=0 


A.14.2. Partes simétrica y antisimétrica 


Como vamos a mostrar a continuación, cualquier tensor se puede descomponer en uno que sea 
simétrico y otro que sea antisimétrico: 
A=As+ÁAa (A.167) 


donde As se refiere a la parte simétrica, 
1 t 
As = 54 +A') (A.168) 


que tiene seis componentes independientes, y Aa se refiere a la parte antisimétrica, 


As=>(4-4') (A.169) 


que tiene tres componentes independientes. 
Sea A un tensor, que expresado en sus componentes simétrica y antisimétrica tendrá la forma: 
1 1 1 1 1 1 
A= A+ A+ Al — A = -(A+4)+-(4-—A! 
2 2 2 2 z. ) 5( ) 
donde el primero y el segundo sumando del último término se corresponden, respectivamente, con un 
tensor simétrico y uno antisimétrico. 


A.15. EJERCICIOS PROPUESTOS 


A.15. Ejercicios propuestos 


Ejercicio 1. Exprese en notación de índices las siguientes cantidades: 


y verifique que AA(BAC)=B(A-C)-C(A-B) 
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Ejercicio 2. Escriba cada una de las expresiones siguientes teniendo en cuenta el convenio de sumación 


de los índices repetidos: 


n dó = Hd + de +- H dl. 
STRE FET EET. 
m ds? = g1 (dx!) + g2 (d? } + g33(dx°}. 


m 1 ad gpadxPdx1. 


Ejercicio 3. Escriba todos los términos de cada una de las siguientes sumas indicadas: 


k 
E a jkX 
qr 
m ApgA 


E Srs 


Ejercicio 4. Halle el ángulo formado por los vectores A = (2,2,—1) y B = (6, —3,2). 
Ejercicio 5. Halle los ángulos que forma el vector A= (3,—6,2) con los ejes coordenados. 


Ejercicio 6. Halle la proyección del vector A = (1, —2, 1) según la dirección de B = (4, —4,7). 


Ejercicio 7. Halle la ecuación del plano perpendicular al vector A= (2,3,6) y que pasa por el extremo 


del vector B= (1,5,3). 


Ejercicio 8. Demuestre que 
|A A B|? + JA- BI? = |A|? |B]? 
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Ejercicio 9. Si A = (3,—1,2), B= (2,1,—1) y Č = (1, —2,2), halle 


(AA B) AC. 


Ejercicio 10. Determine el vector unitario perpendicular al plano formado por A= (2,-6,-3) y 
B= (4,3,—1). 


Ejercicio 11. Siendo A — ail aj + azk y B= bii+ b2j+ bak, demuestre gue 


i j kÈ 
AAB=| aj a aa 
bi b; ba 
Ejercicio 12. Demuestre gue el valor absoluto de 
A- (BAČ) 


es igual al volumen de un paralelepípedo de aristas A, B y C. 
Ejercicio 13. Calcule (27— 37) - [€+ J- $) A(3í—K)]. 


Ejercicio 14. Halle la ecuación del plano formado por los puntos P¡(2,—1,1), P,(3,2,—1) y P3(—1,3,2). 


Ejercicio 15. Siendo ó(x, y,:) = 3x“y — y°z, halle Vo (gradiente) en el punto (1, —2,—1). 


Ejercicio 16. Siendo A = x?zi — 2y°z? j + xy?zk, halle V-A (divergencia) en el punto (1,—1, 1). 


Ejercicio 17. Siendo A =xz%i—2x2yzj+2y2*K halle V AA (rotacional) en el punto (1,—1,1). 


Ejercicio 18. Demuestre que el producto exterior de dos vectores contravariantes A! y Bi, resulta un 
tensor de segundo orden contravariante. 


Ejercicio 19. Sean A;; y A' las componentes de un tensor de segundo orden. Pruebe que A — A; ¡AÑ 
es un escalar. 


Ejercicio 20. Sea Aj un tensor de tercer orden. 


oak ij é š 
= Demuestre que la contracción A;/ es un tensor de primer orden contravariante. 
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= Demuestre que la contracción de los índices i y j produce Ay, el cual no es un tensor. Esto 
muestra que, en general, el proceso de contracción se aplica a índices de distinto nivel. 


Ejercicio 21. Encuentre los autovalores y autovectores asociados a la matriz A: 
1 1 2 
A=| 12 1 
2 1 1 
Muestre que los autovectores son ortogonales. 


Ejercicio 22. Exprese en forma tensorial: 


m la velocidad. 


= la aceleración de una partícula en movimiento. 


Ejercicio 23. Exprese la ley de Newton de la mecánica en forma matricial. 
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